Arbitrarily accurate composite pulse sequences 
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Systematic errors in quantum operations can be the dominating source of imperfection in achieving 
control over quantum systems. This problem, which has been well studied in nuclear magnetic 
resonance, can be addressed by replacing single operations with composite sequences of pulsed 
operations, which cause errors to cancel by symmetry. Remarkably, this can be achieved without 
knowledge of the amount of error e. Independent of the initial state of the system, current techniques 
allow the error to be reduced to O(e^). Here, we extend the composite pulse technique to cancel 
errors to 0(e"), for arbitrary n. 



Precise and complete control over closed quantum sys- 
tems is a long-sought goal in atomic physics, molecular 
chemistry, condensed matter research, with fundamental 
implications for metrology 0, 0] and computation^, 
Achieving this goal will require careful compensation for 
errors of both random and systematic nature. And while 
recent advances in quantum error correction^, 0, Q al- 
low all such errors to be removed in principle, active er- 
ror correction requires expanding the size of the quan- 
tum system, and feedback measurements which may be 
unavailable. Furthermore, in many systems, errors may 
be dominated by those of systematic nature, rather than 
random errors, as when the classical control apparatus 
is miscalibrated or suffers from inhomogeneities over the 
spatial extent of the target quantum system. 

Of course, systematic errors can be reduced simply by 
calibration, but that is often impractical, especially when 
controlling large systems, or when the required control 
error magnitude is smaller than that easily measurable. 
Interestingly, however, systematic errors in controlling 
quantum systems can be compensated without specific 
knowledge of the magnitude of the error. This fact is 
loreQ in the art of NMR, and is achieved using the 
method of composite pulses, in which a single imperfect 
pulse with fractional error e is replaced with a sequence 
of pulses, which reduces the error to 0(e"). 

Composite pulse sequences have been constructed to 
correct for a wide variety of systematic errors 0, 0, 0| . 
These include pulse amplitude, phase, and frequency er- 
rors and can be applied to any system with sufficient 
control. As system control increases, new uses for com- 
posite pulses emerge. A remarkable example is the recent 
teleporation of an atomic state in ion traps [lllIT^ . Bar- 
ret et al. use a composite pulse for individual addressing, 
while Reibe et al. use a composite pulse to perform two 
qubit operations. 

In the context of spectroscopy, the goal is often to max- 
imize the measurable signal from a system which starts 
in a specific state. Thus, while composite sequences have 
been developedjl^J which can reduce errors to 0(e") for 
arbitrary n, these sequences are not general and do not 
apply, for example, to quantum computation, where the 



initial state is arbitrary, and multiple operations must be 
cascaded to obtain desired multi-qubit transformations. 

Only a few composite pulse sequences are known which 
are fully compensating,^^ 0| meaning that they work 
on any initial state and can replace a single pulse without 
further modification of other pulses. As has been theo- 
retically discussed lldll'H.ll?! an d experimentally demon- 
strated in ion traps [ij ll3 and Josephson junctions 
|20j . these sequences can be valuable for precise single 
and multiple-qubit control using gate voltages or laser 
excitation. 

Previously, the best fully compensating composite 
pulse sequence known [ij, [13, 113 could only correct 
errors to 0(e^)[23|. Here, we present a new, and system- 
atic technique for creating composite pulse sequences to 
correct errors to 0(e"), for arbitrary n. The technique 
presented is very general and can be used to correct a 
wide variety of systematic errors. Below, our technique 
is illustrated for the specific case of systematic ampli- 
tude errors, using two approaches. Also discussed is the 
number of pulses required as a function of n. 

The problem of systematic amplitude errors is modeled 
by representing single qubit rotations as 



i?0 (9) = exp 



(1) 



where 6 is the desired rotation angle about about the 
axis that makes the angle (/> with the x-axis and lies in 
the X — y plane, cr^ — cos((/))X + sin(0)y, and X and Y 
are Pauli operators. (9) is the ideal operation, and 
due to errors, the actual operation is, instead, (9) = 
{9{1 + e)), where the angle of rotation differs from 
the desired 9 by the factor 1 + e. Note that cj) and 9 
may be specified arbitrarily, but the error e is fixed for 
all operations, and unknown. 

Two methods for constructing composite pulses. 

A composite pulse sequence i?^' (9) is a sequence of op- 
erations {M^i9)} such that (9) = R^ (9) + 0(e"+i), 

for unknown error e. To construct i?^"' (6*), we begin with 
two simple observations: first, i?^ {—9e) (9) = R^ {9) 
and second, (2fc7r) = ±i?0 (2fc7re) when k is an inte- 
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ger. A composite pulse sequence can thus be obtained 
by finding ways to approximate i?^ (— ^e) by a product 
of operators -R^, (2fc;7re). We obtain this using two ap- 
proaches. 

The first approach we call the Trotter-Suzuki (TS) 
method. Suzuki has developed a set of Trotter formu- 
las that when given a Hamiltonian B and and a series 
of Hamiltonians {Ai} such that B = J2 A-i there exists a 
set of real numbers {pjn} such that 

exp i-iBt) - Yl cxp i-ipjnAit) + 0(r+i) , (2) 

and J2jPjn — 1 _21J. Without loss of generality, we 
may limit ourselves to expansions where the pjn are ra- 
tional numbers, and assume the goal is to approximate 
Ra {~0e). Using Eq. ©, we settle and B = -{9/2)X. 
Then we choose Ai — = rmT{X cos (p + F sin0) and 
A2 = 2r7i7r(X cos — Ysmcj)) where and m fulfill the 
conditions that Amncosip — 6/2 {i.e., Ai+A2 + A:i = B) 
and qjn = Pju^n is an integer. This yields an nth order 
correction sequence 

Fn = nM^(27rgj„)Af_^(47r(7j„)M^(27r(7j„) 
j 

= R„ i-Oe) + 0(e«+i) = 4"' i~6e) (3) 

and the associated nth order composite pulse sequence 
F„Mo (9) = i-ee) Ro {Oe) Rq (0) = i?!,"' (^), thus giv- 
ing a composite pulse sequence of arbitrary accuracy. 

The second approach we refer to as the Solovay-Kitaev 
(SK) method, as it uses elements of the proof of the 
Solovay-Kitaev theorem First, note that rotations 

Uk{A) = I + Ae'' + 0(e*'+^) can be constructed for arbi- 
trary 2x2 Hermitian matrices A, and k > 1, recursively. 
This is done using an observation (from |22||) relating 
the commutator [A, B] = AB — BA to a sequence of op- 
erations, exp(— lAe') exp(— iSe™) exp(i^e') exp(ii?e™) = 
exp([A,B]e'+™)-t-C'(e'+™+i). Thus to generate f/fc (A) it 
suffices to generate C/[fe/2](-B) and U^k/2\{C) such that 
[B,C] ~ A (choices of integers other than [fc/2] and 
[fc/2j which sum to k are also fine, but less optimal). 

Next, we inductively construct a composite pulse se- 
quence Fn for Rq {9). Note that the first order correction 
sequence can be written as Fi — (27r) M-^ (27r) = 
Ro{~9e) + 0{e^) by selecting 47rcos((/)) = 9. Assume 
we have F„ = Rq {-9t) - iA„+ie"+i + 0(e"+2). We 
can then construct a sequence to correct for the next or- 
der, using F„+i = [/„+i(A„+i)F„, where L/„+i(A„+i) is 
constructed as above. Iteratively applying this method 
for k = 1 , . . . , n yields an nth order composite pulse se- 
quence, FnMo{9) = i?!)"' (6*), for any n. This method, 
which appears to be unrelated to previous composite 
pulse techniques H, 0|, gives an efficient algorithm to 
calculate sequences for specific 9 and (/) but not necessar- 
ily a short analytical description of the sequence. Fur- 
thermore, the Solovay-Kitaev technique relies on general 



properties of Hamiltonians and can be applied without 
modification to other systematic error models, e.g., fre- 
quency errors. 

Examples. The TS and SK techniques described above 
are general and apply to a wide variety of errors; ex- 
plicit application of the techniques to generate Rq^^ (9) 
sequences for specific n can take advantage of symmetry 
arguments, composition of techniques, and relax some of 
our assumptions to minimize both the residual error and 
the sequence length. 

First, we explicitly write out the TS composite pulses 
and connect them to the well-known pulse sequences 
of Wimperis We choose to use the TS formulas 

that are symmetric under reversal of pulses, i.e an 
anagram. These formulas remove all even-ordered errors 
by symmetry, and thus yield only even-order composite 
pulse sequences. For convenience, we introduce the no- 
tation 5i((/)i, 02, — Af^j (m7r) (2m7r) (tott) 
and Sn{4>i,4>2,m) = S'„_i(0i, 02, to)''" ' x 
5„_i(0i, 02, -2m)5„_i(0i, 02,to)''" . We can now 
define a series of n order composite pulses Pn as 

PO = Mo {9) (4) 
P2 = M^, {27t) M_^, (47r) M^, (2^) PO (5) 
P2j = 5,(0„-0„2)PO (6) 

where 0j = cos^^ ^'^^ ~ (S^^-*^^^ — 2)/j_i when 

/i = 1. P2 is exactly the passband sequence PBl de- 
scribed by Wimperis 0. Fig. ^ compares the perfor- 
mance of these high-order passband pulse sequences. 

Wimperis also proposes a similar broadband sequence, 
BBl=5i(0Bi,30Bi,l)PO where 0si = cos"! (-|^). 
The broadband sequence corrects over a wider range of 
e by minimizing the first order commutator and thus the 
leading order errors. Furthermore, although BBl and 
PBl appear different when written as imperfect rota- 
tions, a transformation to true rotations shows that they 
have the same form, 

PBl = M^,{2tt)M^^,{A^)M^,{2tt)PQ 

= R4,, (2^e) (47re) R^, (2^e) PO 

BBl = M^,, (^) M3^,, (2^) M^,, (tt) PO (7) 

= R^,, (Tre) (2^e) R^,, (^e) PO. (8) 

This "toggled" frame suggests a way to create higher- 
order broadband pulses. One simply takes a higher- 
order passband sequence and replaces each element 
Si{(j)j,~(f)j,m) with Si{4>Bj,-4'Bj + A(t)Bj{rn/2 
mod 2), to/2) where satisfies the condition 

cos(0_Bj) — 2cos(0j). Applying this to Pn creates 
a family of broadband composite pulses, Bn. 

Similar extensions allow creation of another kind of 
composite pulse (useful, for example, in magnetic reso- 
nance imaging), which increase error so as to perform 
the desired operation for only a small window of error. 
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FIG. 1: Comparison of the narrowband and broadband com- 
posite pulse sequences generated by the TS method. The 
Wimperis BBl, PBl, and NBl sequences are included in this 
family, and are equivalent to B2, P2, and N2. 



Such "narrowband" pulse sequences Nn may be obtained 
starting with a passband sequence, Pn, and dividing the 
angles of the corrective pulses by 2. These higher-order 
narrowband pulses may be compared with the Wimperis 
sequence NBl 15], as shown in Fig.^ 

The SK method yields a third set of nth order com- 
posite pulses, SKn, and for concreteness, we present an 
explicit formulation of this method. It is convenient to 
let f7„x(a) = / - ia"f e" -I- 0(e"+i), such that one can 
then generate J7„z(o) = Afgo (-7r/2) C/„x(a)M9o (7r/2) 
and [/^^(a) = M45 (7r/2) [/„x (0)^45 (-7r/2) J^]. Using 
the first-order rotations 



a 



M-, 



(2n 



a 

in 



(9) 



where cj) — cos ^(a/(47r [^] )), as described 
above, we may recursively construct Unx{a) — 

U\n/2\ Y {a)U \n/2^ Z («) U L„/2J Y (-«) C/p„/2] Z (-fl) , for 

any n > 1 and any a. 

With these definitions, the first order SK composite 
pulse for i?!)"' [9) is simply 

SKI = Uix{e)Mo (e) = Ro {6) - i^e^ + 0{e^) . (10) 

From the 2x2 matrix A2, we can then calculate the 
norm IIA2II and the planar rotation that performs 
i?A2 (~^2)^]4^ — 11^2 11-^- The second order SK compos- 
ite pulse is then 



SK2 = MX'U2x{\\A2\\'^^)Ma,SK1 



Ro (9) - z^e^ + 0{e*) 



(11) 
(12) 



where is the imperfect rotation corresponding to the 
perfect rotation Ra2 ■ 



The nth order SK composite pulse family is thus 

SKn = M^^i[/„x(P„ir/")MA„SK(n-l) (13) 
^ i?o(0)_j^^e»+i+O(e"+2). (14) 

A nice feature of the SK method is that when given 
a composite pulse of order n described by any method 
one can compose a pulse of order n 4- 1. The "pure" 
SK method SKn is outperformed in terms of both error 
reduction and pulse number by the TS method Bn for 
n < 4. Therefore, we apply the SK method for orders 
n > 4 using B4 as our base composite pulse. We label 
these pulses SBn. 

Performance and efficiency. Two important issues 
with composite pulses are the actual amount of error re- 
duction as a function of pulse error, and the time required 
to achieve a desired amount of error reduction. These 
performance metrics are shown in Fig. |21 comparing the 
SK, broadband, and passband composite pulses for vary- 
ing error e, and = 0, and using as the composite pulse 
error E = \\R^ (6) - i?^"' (6*) ||. We find that for prac- 
tical values of error reduction, n < 30, the number of tt 



pulses required to reduce error to 0{e^' 
which is close to the lower bound of ' 



grows as - 



-,3.09 



which can be 



analytically derived|23|. In contrast, the TS sequence Bn 
requires 0(exp(n^)) pulses. 

For a wide range of base errors e, the TS formulation 
out-performs the SK method in achieving a low com- 
posite pulse error, E. The recursive nature of the TS 
methods builds off elements that remove lower order er- 
rors, resulting in a rapid increase of pulse number and a 
monotonic decrease in effective error at every order for 
any value of the base error. However, the SK approach 
is superior to the TS method for applications requiring 
incredible precision, E < 10^^^, from relatively precise 
controls, e < 10~^. 

The SK and TS pulse sequences presented here are 
conceptually simple but may not be optimal. Integrating 
ideas from both methods, we can develop new families 
of composite pulses. As an example, the SK method re- 
lies on cancellation of error order by order by building 
up sequences of 27r pulses. However, there is no reason 
that the basic unit should be a single pulse. Instead, one 
can build a sequence from TS (B2) style pulse triplets, 
G{(j)i) = Si{(f>i,34>i,l). By using an additional symme- 
try that the Tr (rG(-(^i)G'((/)i)) = 0, the leading order 
error is guaranteed to be proportional to X at the cost 
of doubling the pulse sequence. The resulting pulses are 
of length exp(n) (compared to exp(n^) for TS), broad- 
band compared to SK sequences, and described in detail 
in "2^. 

Conclusions. We have presented a set of tools that al- 
lows one to generate arbitrarily accurate composite pulse 
sequences for systematic, but unknown, error. As an 
example, we have constructed explicit composite pulse 
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FIG. 2: Composite pulse error, E, as a function of base er- 
ror, e for a variety of composite pulse sequences. Pn, Bn, 
SKn and SBn are the nth order passband, broadband, SK, 
and combined B4-SK sequences, respectively. The number in 
the brackets refers to the number of imperfect 2ti rotations in 
the correction sequence. Note how pulses of the same order 
(such as P6, B6, SK6, SB6) have the same slope (asymptotic 
scaling) for low values of e, but can have widely varying per- 
formance when e is large. The inset plots the scaling of this 
sequence length with order n for SKn (SBn is very similar) 
for n < 30 and compares it with the upper bound obtained 
with numerical methods. 

sequences for errors in rotation angle. These can be 
constructed with 0{rt') pulses, for n < 30. For high- 
precision applications such as quantum computation, 
these pulses allow one to perform accurate operations 
even with large errors. Practically, the B4 and B2=BB1 
pulse sequences seem most useful, depending on the mag- 
nitude of error. 

While we have focused on composite pulse sequences 
for rotation errors, we emphasize that these meth- 
ods also apply to correcting systematic errors in con- 
trol phase and frequency [23. For example, a fre- 
quency error can be represented for an expected rota- 
tion Ro{e) as M^{e) = exp {-i{e/2X + \e/2\5Z)) . Note 
that M^(6»/2)M^(-6')M^(6i/2) yields to first order in 5 the 
phase shift, Uiz{29S). Starting with any fully compen- 
sating composite pulse sequence that corrects frequency 
errors to order 5^, e.g. CORPSE^, and the basic op- 
eration Uiz{29S), one can then apply the SK technique 
to create a pulse sequence of 0(^"l|2^. 

Furthermore, the TS and SK approaches can be ex- 
tended to any set of operations that has a subgroup iso- 
morphic to rotations of a spin. For example, Jones has 
used this isomorphism to create reliable two qubit gates 
based on an Ising interaction to accuracy 0{e^) Sim- 
ilarly, the techniques outlined here can immediately be 
applied to gain arbitrary accuracy multi-qubit gates. In- 



terestingly, the TS formula can be directly applied to any 
set of operations, if the operations suffer from propor- 
tional systematic timing errors. Therefore, this control 
method could also be applied to classical systems. 
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